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McCOIVAN’S  SOLITARY  WAVE  EXPANSIONS 


ABSTRACT 

In  1894  McCowan  proposed  that  the  structure  of  solitary  waves  be 
represented  by  a  certain  series  expansion;  but  he  gave  an  estimate  for 
only  the  first  term  of  this  expansion.  Any  truncation  of  such  an  expan¬ 
sion  compromises  the  essential  ideal  of  the  steadiness  of  a  solitary 
wave;  also,  the  approximations  which  he  adopted  might  not  be  the  best 
ones  for  the  diverse  particular  purposes  for  v/hich  this  wave  may  be 
used.  In  this  paper  a  more  precise  structure  of  this  wave  is  constructed 
by  implementing  McGowan's  expansion  to  higher  orders.  The  coefficients 
in  this  expansion  are  determined  by  minimizing  a  certain  reasonable 
measure  of  the  imperfection  of  this  expansion.  This  minimization  is 
done  by  using  the  ENMIN  subroutine. 
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INTRODUCTION 

In  the  Theory  of  Water  Waves  there  often  arises  the  concept  of  a 
solitary  wave,  often  spoken  of  as  a  *  soli ton'.  It  is  a  single,  two- 
dimensional  intumescence  that  propagates  (over  a  level  bottom,  and  with 
a  series  of  traditonal  simplifying  assumptions)  supposedly  without  any 
change  in  its  profile  and  velocity;  its  velocity  depends  upon  its  in¬ 
tensity.  No  exact  and  concise  expression  for  the  structure  of  such  a 
wave  is  known.  Until  recently,  all  theories  of  this  wave  have  been  in 
the  form  of  expansions;  if  these  theories  may  have  appeared  to  be  simple, 
this  is  so  only  because  for  various  particular  uses  of  this  concept 
(and  specifically,  for  solitary  waves  of  low  intensity)  the  first  terms 
of  such  expansions  have  been  considered  sufficient.  Any  truncation  of 
such  expansions,  of  course,  possesses  shortcomings.  Thus,  the  popular 
Boussinesq  model  of  solitary  waves  compromises  (as  shown  by  Munk  ) 
the  generally  presumed,  Laplacian,  governing  equation;  while  the  com¬ 
monly  quoted  first-order  variant  of  McCowan's  theory  compromises  the 
dynamic  boundary  condition  on  the  pressure  at  the  surface.  There  have 

been  a  number  of  attempts  to  reduce  such  shortcomings  by  resorting  to 
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different  expansions  ,  or  to  higher-order  terms  of  such  expansions 
(with  the  coefficients  of  such  terms  determined  analytically);  but  all 
such  attempts  drastically  complicate  these  theories,  Byat-Smith4  and 
Strelkofr  have  recently  proposed  a  radical  resolution  of  such  short¬ 
comings  by  obtaining  the  profile  of  this  wave  numerically;  but  this  ap¬ 
proach  involves  even  more  sophisticated  and  laborious  procedures  (in¬ 
volving  the  solution  of  an  integro-differential  equation),  is  not  im¬ 
mune  to  'numerical1  (truncation-and-roundoff)  errors,  and  requires  much 
further  work  if  the  interior  flow  field  is  desired. 

For  certain  practical  purposes  such  as  those  of  Reference  6,  where 
in  the  representation  of  a  solitary  wave  it  is  desired  to  select  a  good 
compromise  between  simplicity  and  precision  (or  at  least  to  inspect, 
and  to  exhibit,  the  effectiveness  of  various  tradeoffs  between  these 
desiderata)  the  most  efficient  course  appeared  to  be  that  of  extending 

References  are  found  on  p-ye  21. 
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McCowan's  theory  to  higher  orders,  with  the  coefficients  of  the  higher- 
order  terms  determined  ’empirically'  by  computations. 

GENERAL  EXPANSION  PROPOSED  BY  McCOWAN 

McCowan  originally  proposed  that  the  complex  stream  function  of  a 
solitary  wave  (in  a  stationary  coordinate  system,  at  the  instant  when 
the  o’est  of  the  wave  passes  the  plane  of  x  =  0)  be  represented  as 

♦  +  =  I  a.  tan-*  [~(z  +  ix)],  j  =  1,  3,  5...  (1) 

where  x  and  z,  respectively,  are  the  horizontal  and  vertical  distances 
to  a  water  particle  from  the  point  at  the  bottom  under  the  wave's  crest; 
the  real  stream  function  < P  and  the  real  velocity  potential  <{>  are  here 
defined  by  the  fact  that  the  horizontal  and  vertical  components  of  the 
water  velocity  are 


u  s  a$/*x  =  atydz,  v  =  dt'/oz  =  -  a^/ax,  (2) 


j"d  the  real  constants  and  m,  as  well  as  the  wave's  celerity  C,  re¬ 
main  to  be  determined.  Incidentally,  the  parameter  m  can  to  some  ex¬ 
tent  be  specified  arbitrarily,  since  it  indicates  the  intensity  of  the 
wave  (m  is  in  the  nature  of  a  reciprocal  of  some  "effective’1  length  cf 
the  wave,  and  the  more  intense  wave  is  the  'more  concentrated',  or  the 
shorter,  one).  The  intensity  of  the  wave,  however,  ma>  more  conven¬ 
iently  be  specified  by  the  ratio  y  of  waveheight  to  undisturbed  depth; 
thus  a^ ,  m,  and  C  ought  to  be  expressed  in  terms  of  y.  It  might  be  no¬ 
ticed  that  the  notation  of  (1)  and  (2),  particularly  well  suited  to  this 
problem,  differs  slightly  from  what  would  be  more  customary  today. 


McCowan,  however,  at  once  in  effect  confined  his  theory  to  the 
first  term  of  this  expansion  -  by  using  a  couple  of  ingenious  tricks 
for  the  choice  of  the  coefficients  a^ ,  m  and  C.  The  difficulties  which 
he  faced  in  his  ground-breaking  1894  proposal  were  precisely  of  the  type 
that  are  readily  resolvable  with  today's  computing  technology. 
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REQUIREMENTS  UPON'  Till;  COEFFICIENTS  IN  THE  EXPANSION 

An  obvious  reason  for  McCowan's  choice  of  the  basic  function  (the 
tangent  of  the  complex  coordinate  in  a  certain  scale)  in  this  expansion 
is  that  the  first-order  term  alone  already  yields  a  configuration  that 
indeed  is  very  similar  to  the  ideal  described  by  Russell  empirically, 
and  had  in  mind  by  Boussinesq  and  Rayleigh.  The  reason  for  McCowan's 
providing  -  at  least  in  principle  -  the  infinite  number  of  coefficients 
3^  can  be  easily  seen  if  we  consider  the  requirements  which  a  solitary 
wave  must  satisfy. 

The  resort  to  an  analytic  function  of  a  complex  coordinate  at  once 
satisfies  thr  governing  (Laplace's)  equation.  The  choice  of  a  function 
(the  tangent)  such  that  its  real  part  vanishes  for  z  =  0  satisfies  the 
bottom  boundary  condition.  Fhc  remaining  conditions  are  best  inspected 
ip.  a  coordinate  system  moving  with  the  wave,  in  which  the  complex  stream 
function  becomes,  say 

T  +  i4>  =  (i^+i <J>)  -  C(z+ix)  (3) 

so  that  transfer  to  this  system  reduces  u  by  C  and  leaves  v  unchanged. 

In  particular,  the  kinematic  boundary  condition  on  a  wave's  surface  (the 
impenetrability  of  the  free  surface)  then  can  be  satisfied  simply  by  defi¬ 
nition:  a  streamline  can  be  taken  as  the  profile  of  the  wave.  The  dy¬ 
namic  boundary  condition  (the  continuity  of  pressure  across  the  air-water 
interface)  is  further  simplified  because  McCowan's  theory  is  one  of  those 
which  neglect  the  density,  and  hence  the  pressure,  of  air.  This  condi¬ 
tion  then  reduces  to  p  =  0,  and  Bernoulli's  equation  for  a  water  particle 
at  the  surface  takes  the  form 

2 

«I>  +  Q  /2  +  gZ  =  constant 
2 

=  C  /2  +  gh  as  x  *  ±  ~ 
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where  h  is  the  undisturbed  depth,  Q“  is  (U  -  C)  +  V  ,  nnd  Z,  U,  V  are 
the  surface  values  of  z,  u,  v.  For  any  instantaneous  functions  Z(x) 
and  vj(x)  this  equation  merely  yields  4>t;  but  the  solitary  wave  is  dis¬ 
tinguished  from  other  waves  by  still  another,  the  fifth,  condition,  that 
of  its  steadiness:  viz.,  by  the  fact  that  its  4*  must  vanish  identi¬ 
cally  everywhere.  It  is  in  order  to  assure  such  vanishing  tnat  we  have 
been  given  the  freedom  of  choosing  m,  C  and  the  a^'s. 

A  CONVENIENT  CHOICE  OF  UNITS  FOR  ANALYSIS 

If  h  and  the  unknown  C  are  •  iken  as  units  of  distance  and  velocity 

2 

(so  that  the  unit  of  time  is  h/C;  of  acceleration,  C  /h;  and  of  S',  Ch) , 
this  C  drops  out  of  our  analysis,  but  can  easily  be  recovered.  Tn  lieu 
of  C  there  appears  another  unknown,  the  acceleration  of  gravity  in  these 
units,  viz.,  the  quantity  gM  =  g/(C2/h);  but  this  gM  happens  to  enter 
the  theory  in  such  a  tanner  that  it  is  much  more  easily  determinable  than 
C.  Another  interpretation  of  gM  is  convenient:  writing  its  definition 
as 


R  2  (gV1/2  =  C/(gh)1/2  (5) 

1/2 

and  recollecting  that  (,h)  is  the  celerity  of  Airy  low-and-long  (or 
"shallow-water")  waves  of  small  amplitude,  (5)  is  seen  to  correspond  to 
the  definition  of  Mach  number  in  Compressible  Fluid  Mechanics  -  whereby 
the  solitary  wave  corresponds  to  a  shock  wave.  The  term  "Russell  num¬ 
ber",  after  the  discoverer  of  this  long-unrecognized  wave,  is  suggested 

for  R.  It  should  be  noted  that  in  both  the  1844  empirical  concepts  of 

Scott  Russell,  and  in  the  1871  Boussinesq  theory  it  turns  out  that 

R  -  (1+v) 1/2  «>) 
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-  a  v;lue  which  our  extension  of  McCowan's  theory  modifies  only  very 
slightly  (while  in  the  original  tricky  version  of  his  first-order  theory 
R  appeared  to  be  much  smaller). 

The  quantities  m,  a^  and  a^  for  j  =  3,  S...  in  these  units  will  be 
denoted  by  M,  N/M,  and  AuN/M,  where  N  turns  out  to  be  a  convenient  fac¬ 
tor  for  all  velocities.  If  the  tangent  in  (1)  is  T,  the  complex  stream 

2  4 

function  appears  as  (NT/M)  (1+A_T  +A.T  +...).  Our  problem  thus  reduces 

3  5  M 

to  a  determination  of  such  M,  N,  A_,  Ac...  and  g  as  would  best  anni- 
hi  late  4>t. 

CRITERION  FOR  THE  IMPERFECTION  OF  A  SOLITARY  WAVE  MODEL 

Given  a  tentative  set  of  parameters  M,  N,  g^1,  A^,  A^...  of  a  soli¬ 
tary  wave,  we  may  compute  fin  a  manner  that  will  presently  be  shown) 
for  every  x  the  quantity 

D  =  2+t  =  1  -  Q2  -  2g'M(Z-l)  (7) 

which  we  shall  consider  as  a  discrepancy,  or  unsteadiness.  Clearly,  in 
order  to  formulate  a  criterion  that  would  characterize  this  tentative 
set  of  parameters  without  specifying  x  we  need  merely  to  compute  D  for 
several  values  of  x  and  do  some  averaging.  The  point  is  that  the  man¬ 
ner  of  such  averaging  is  essentially  arbitrary.  In  particular,  we  find 
it  both  informative  and  practical  to  consider,  and  to  minimize,  the 
positive-definite  quantity 

rY  2 

F  =  /  D  a  Y  (8) 

•'0 

where  Y  =  Z  -  1.  By  viewing  D  as  D(Y),  and  thus  taking  a  set  of  points 
uniformly  spaced  in  Y,  we  give  we; t  mainly  to  points  in  the  main  body 
of  the  wave,  and  give  only  small  wc!„nr  to  the  points  in  the  "far 
tails";  also,  we  take  advantage  of  the  fact  that  at  Y  =  0  (viz.,  at 
x  =  ±  »)  D  vanishes.  In  this  connection  it  should  be  noted  that  one 
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of  McGowan’s  ’tricks’  was  the  basically-arbitrary,  and  unnecessarily 
stringent,  requirement  as  to  the  manner  in  which  I)  -<•  0  as  x  >•  •<  <»:  he 
required  that  the  curve  of  D(Y)  be  tangent  to  the  Y-axis  at  the  origin 
(at  Y  -v  0).  He  used  this  condition  to  determine  the  celerity  C  of  trio 
solitary  wave.  Unfortunately,  this  method  amounted  to  determining  C 
from  the  behavior  of  the  "far  tails"  (which  may  be  of  little  interest  in 
such  a  context  as  that  of  Reference  6);  and  to  disregarding  the  rela¬ 
tively  large  discrepancies  D  which  thereby  result  in  the  main  body  of 
the  wave.  If  his  requirement  on  the  manner  of  the  vanishing  of  D  at 
Y  -  0  is  relaxed,  much  smaller  P  can  be  achieved,  even  when  only  the 
first  term  of  (1)  is  used;  and  in  particular,  this  reduction  of  P  means 
that  the  dynamic  boundary  condition  is  more  nearly  satisfied  in  the  main 
body  of  the  wave. 

COMPUTATION  OP  Till-  lMPURPHCTION  CiUTP.IUON 

Computation  of  D  for  use  in  (8)  requires  a  determination  of  x  for 
a  prescribed  Z,  vis.,  the  determination  of  the  wave's  profile;  for  once 
x(Z)  is  known,  U,  V  and  Q  can  be  gotten  by  use  of  (2).  In  the  moving 
coordinate  system  the  ieal  stream  function  is  Y  -  -  z;  and  at  the  sur¬ 
face,  at  x  4  1  “  (where  z  ~  l  -*■  1  and  ij»  *>  0)  this  Y  becomes  -  1.  Wave's 

profile  thus  i?  the  streamline  of  Y  =  -1,  whereby  i{i=Y  +  Z=  Z-  1  =  Y, 

Thus  by  equating  the  real  part  of  (1)  to  Z  -  1,  viz.,  by  writing 

*  (x,Z)  =  Y  (9) 

where  Z  and  Y  are  prescribed,  we  have  an  equation  for  determining  x(Y). 

Th;  two  end  points  in  the  range  0  £  Y  £  y  afford  simplifications;  x(y1 

is  simply  0,  while  x(0)  is  of  no  interest  since  it  yields  D  =  0.  For 
other  values  of  Y  the  equation  (9)  can  readily  be  solved  by  the  Newton¬ 
ian  iterations.  In  such  iterations  the  first  approximation  for  x(Zl  can 
of  course  be  taken  as  that  obtained  with  seme  preceding  tentative  set  of 
the  parameters  M,  N,  gW,  A^,  Aj. . . . ;  incidentally,  V  appears  as  a  byproduct 
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ot  these  iterations*  With  x(Z)  known,  the  determination  of  U  by  (2),  D 
by  (7)  and  F  by  (8)  is  straightforward. 

Hie  minimisation  of  F  -  and  hence,  the  determination  of  the  best 
set  of  the  wave’s  parameters  for  a  given  y  -  is  in  practice  made  possi¬ 
ble  by  the  existence  of  the  FORTRAN  subroutine  FNM1N  (’function  minimi¬ 
zation  without  taking  the  derivatives’).  The  value  of  F  computed  for  a 
tentative  set  of  parameters  is  simply  returned  to  that  subroutine  - 
which  then  decides  what  alternative  combination  of  parameters  is  to  be 
tried  next,  and  the  search  continues  until  the  subroutine  is  satisfied, 
viz.,  until  all  parameters  are  optimized  to  a  prescribed  precision. 

STRUCTURE  OF  THE  REFINED  SOLITARY  WAVE 

Results  of  our  minimization  of  F  can  be  expressed  by  a  statement 
of  the  values  of  parameters  M,  N/M,  A^,  A,,...  and  R  for  a  series  of 
values  of  y.  Figure  1  shows  the  behavior  of  the  function  D(Y)  for  dif¬ 
ferent  models  of  a  wave  of  y  =  .5.  Table  I  shows  the  wave  parameters 
for  this  y  for  a  series  of  models  of  this  wave,  and  also  shows  the  pre¬ 
cision  criterion  F,  and  the  number,  n,  of  intervals  into  which  the  range 
of  Y  has  been  divided  in  this  computation.  The  expansion  (1)  is  seen 
to  be  not  of  the  ’orthogonal’  type,  in  that  the  passage  to  the  next- 
higher-order  model  does  modify  the  coefficients  that  enter  into  the  pre¬ 
ceding  model;  but  there  is  a  modicum,  so  to  say,  of  ideal  orthogonality 
in  that  such  modification  is  relatively  minor.  The  inclusion  of  a  jth 
term  in  (i)  seems  to  contribute  a  ripple  of  (j+l)/2  humps  that  attempts, 
so  to  say,  to  fill  in  the  imperfections  in  the  profile  of  the  wave  that 
are  left  by  the  (j-2)  th-order  model.  Table  II  shows  the  coefficients 
in  (1)  for  the  1st,  3rd,  and  5th  order  models.  An  inspection  of  this 
table^  suggests  that  for  the  larger  values  of  y  a  greater  number  of 
terms  in  (1)  ought  to  be  considered;  in  particular,  it  has  shown  that 
for  y  =  .5  the  expansion  (1)  ought  to  be  extended  at  least  to  the  9th 
order.  Table  III,  in  particular,  shows  that  the  celerity  of  the  soli¬ 
tary  wave  is  much  closer  to  its  classical  estimate  by  (6)  than  to 
McCowan’s  estimate. 
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The  background,  the  theory,  the  limitations  -  and  further  numerical 
illustrations  -  of  this  extension  of  McGowan's  theory  of  solitary  waves 
are  given  in  greater  detail  in  Reference  6.  In  particular,  we  purposely 
refrain  here  from  considering  the  much-discussed  question  of  the  esti¬ 
mate  of  the  upper  limit  of  y  for  the  solitary  wave.  For  the  higher  y 
the  crest  of  this  wave  sharpens,  and  the  approach  along  the  lines  of 
(1)  becomes  awkward;  but  there  are  other  reasons  for  believing  that  such 
a  limit  is  a  phantom,  for  most  of  the  classical  assumptions  postulated 
in  all  theories  of  solitary  waves  would  cease  to  be  applicable  long  be¬ 
fore  such  a  limit  is  approached. 
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(  For  Tails  )  (Crest) 


Figure  1.  The  variation  of  the  discrepancy  D,  equation  (7), 
for  different  models  of  solitary  wave  of  height 
y  ~  .  5 . 
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